Semiclassical theory of the emission properties of wave-chaotic resonant cavities 
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We develop a perturbation theory for the lifetime and emis- 
sion intensity for isolated resonances in asymmetric resonant 
cavities. The inverse lifetime P and the emission intensity 
1(9) in the open system are expressed in terms of matrix el- 
ements of operators evaluated with eigenmodes of the closed 
resonator. These matrix elements are calculated in a semi- 
classical approximation which allows us to represent V and 
1(9) as sums over the contributions of rays which escape the 
resonator by refraction. 

PACS numbers: 42.55.Sa, 05.45. Mt, 42.25.-p 

The semiclassical theory of wave equations with non- 
integrable classical (short wavelength) limits has been 
a topic of great interest over the past decade. Much 
progress has been made for closed systems based on pe- 
riodic orbit theory, particularly in the case of fully chaotic 
systems with no stable periodic orbits or families of reg- 
ular orbits (KAM tori) M . One may for example obtain 
eigenvalues with an energy resolution inversely propor- 
tional to the length of the orbits summed [Q. Much cur- 
rent work addresses the generic case of quantum/ wave 
systems with mixed dynamics in the classical limit, i.e. 
with stable, unstable and regular regions in phase space 
H . Here we focus on the situation where a system with 
mixed dynamics is coupled to an infinite asymptotic re- 
gion and the eigenstates become resonances with complex 
eigenvalues A crucial quantity to calculate in this 

case is the resonance width (imaginary part of the eigen- 
value) . A basic problem here for semiclassical approaches 
is that the resonance width may be exponentially smaller 
than the level spacing, making a direct periodic orbit 
method very difficult. 

Although our approach could be applied to quantum 
hamiltonians, we will focus here on an important ap- 
plication for such a theory: the formally similar prob- 
lem of dielectric microcavities. Such resonators are used 
as high-Q resonators for microlasers and optical spec- 
troscopy [^U^]. If such resonators are rotationally sym- 
metric the wave equation is separable and analytic so- 
lutions for resonance properties are possible; moreover 
the corresponding ray dynamics (within the resonator) 
is integrable with conserved angular momentum. When 
such resonators are deformed from rotational symmetry 
the ray dynamics is generically mixed and no analytic 
results were known for resonance properties, except for 
very small deformations which can then be treated per- 
turbatively ||. However some time ago J§-[l^] it was 



noted that such asymmetric resonant cavities (ARCs) are 
equivalent to quantum billiards with refractive escape 
and a ray-optics model was developed for the lifetime 
and emission pattern. Although quite successful for low- 
index resonators where only whispering gallery modes 
are important, this ray model has several limitations. 1) 
It only describes quantitatively the whispering gallery 
modes of such resonators. 2) It neglects all interference 
effects. 3) It is valid only near a certain adiabatic limit 
for the phase space flow |jLljT^| . However recent experi- 
mental results on deformed cylindrical semiconductor mi- 
crolasers (a high-index, high-gain system) demonstrated 
high-power and highly directional emission from "bow- 
tie" modes which are not well described by this model 
Jl3] | . Therefore there is both general theoretical and spe- 
cific experimental motivation for the semiclassical theory 
of ARCs which we present below. 

The key simplification we employ to make the prob- 
lem tractable is to assume that the eigenstates of the 
"closed" system are known and develop a theory which 
is perturbative in the coupling to the outside but not 
in the deformation. The resulting analytic formulae are 
convenient to use in combination with either numerical 
or semiclassical calculations for the closed problem, and 
have a simple physical interpretation in terms of refrac- 
tively escaping rays. For the case of modes associated 
with stable orbits of the closed system a full analytic 
semiclassical solution can be found which depends only 
on nkRo (k is the wavenumber of the radiation, Rq is the 
average radius of the ARC and n is the index of refrac- 
tion) and properties of the classical orbit (e.g. its period 
and stability). Finally we note that the numerical solu- 
tion of such problems becomes intractable due to either 
instabilities or convergence problems fll4| when kRo is of 
order 100, well below the range of many optical exper- 
iments. In this case the semiclassical theory presented 
here is (to our knowledge) the only available method of 
solution. 

We consider a long cylindrical dielectric resonator 
of uniform index n with a non-circular cross-section 
parametrized by the distance R(4>) between the origin 
and a point on the the boundary in the ^-direction. 
We consider the TM polarization (assuming k z — 0) for 
which the electric field is perpendicular to the cylinder 
axis and both the field and its derivative are continu- 
ous at the boundary. A quasibound state with complex 
frequency u> = nek can be written in the form: 
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where 



a m H+ (nkr) + j3 m H m (nkr) if r < R(4>), 
jmH^ (kr) otherwise. 



where are Hankel functions. Such resonant solutions 
can only exist for complex k as there is no incoming wave 
for r > R{4>); the imaginary part T = — Imk is the reso- 
nance width. 

Defining coefficient vectors |a) , |/3) , I7) from Eq. ([!]), 
the regularity of the wavefunction at r = requires \a) — 
. Using the continuity conditions to eliminate I7) and 
the resonance condition can be expressed as the non- 
hermitian eigenvalue problem 



[H (k) + V(k)]\a) = 0, 



(2) 



where 7io and V are the hermitian and antihermi- 
tian part, respectively, of the matrix Ai = (J' J') — 
(l/n)(J'H+') (H-H+y 1 (H-J). Here, M is defined in 
terms of matrices (ZZ) with matrix elements 



{ZZ) mi = — I d<P Z m (k )Z e (k) e 



(3) 



where Z m (k) and Z rn (k) stand for either (kR(<fi)), 
the Bessel function J m (nkR (</>)), or their derivative, and 
k = fco — iT. 

The perturbation theory is based on the observation 
that for a narrow resonance r <C A, (where A is the 
resonance spacing), the antihermitian term V is a small 
correction in Eq. (||) which may be expanded in powers 
of r, to obtain to leading order in V: 



\V\aW) 



(a(°) \{dHo/dk)\a^)y 



(4) 



where the matrix elements are taken at the real wave 
number fc defined by the hermitian eigenvalue problem 



Ho (fco) |a (0) ) = 0. 



(5) 



The operator Hq describes a closed resonator with 
boundary conditions to be specified below. Note, that 
the denominator in Eq. (^) depends only on the prop- 
erties of the closed system, and can be regarded as a 
normalization factor. 

Similarly to (jj), one can express the far-field emission 
intensity from the resonance in terms of the eigenstates 
of H , 



no) = 1 E z u ' e 



(6) 



These results represent a general perturbation theory 
for the problem requiring only small T and may be of 



value in their current form. However, far more interest- 
ing results can be obtained by taking advantage of the 
semiclassical limit in which the wavelength of the radia- 
tion is assumed much smaller than the average radius of 
the resonator nkRo 3> 1. This condition is satisfied in 
essentially all optical experiments. In this limit Hankel 
and Bessel functions may be well-represented by the " ap- 
proximation by tangents" jl5| and the matrix elements 
and angular momentum sums in (Q), (^|) can be evalu- 
ated in the stationary phase approximation. In this ap- 
proximation, the boundary condition corresponding to 
(El) becomes 



i S ' fc„ 



exp (iip) \m)(m\ 



\E) = 0, 



(7) 



where \m) is an eigenfunction of angular momentum 
about the cylinder axis, V„ is the normal derivative at 
the boundary of the resonator and k n (m,(f>) the normal 
component of the wavevector k, and (p(k n ) the Fresnel 
phase change upon reflection from a flat dielectric bound- 
ary. When n ^> f, Eq. (Q) reduces simply to von Neu- 
mann boundary conditions. 

We first present the results for 1(6). Evaluating Eq. 
(^) in the stationary phase approximation gives three 
conditions: 1) The matrix (H~ H + )~ 1 is diagonal. 2) 
An internal angular momentum m only couples to the 
angular momentum I which allows the corresponding ray 
to satisfy Snell's law (either in transmission or reflection). 
3) This ray with angular momentum £q must be emitted 
into angle 9 in the far-field (see Fig. la). The last two 
conditions imply that there exists at most one such ray 
with angular momentum to for each value of to, 9, kR. If 
to corresponds to an angle of incidence sin^ ~ raj nkR 
which is totally-internally reflected (m/kR > 1) then 
there is no such ray, and the corresponding component 
of the field a m does not contribute to emission, hence 
1(9) is only due to refracted rays. At large index the 
contribution from evanescent escape (tunneling) will be 
important and is not described by this level of approx- 
imation. One finds for the (non-normalized) emission 
intensity 
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I (9) = £ m exp(im0±) e m (6) a 



(0) 



(8) 



where e m (9) = J2± B ( ( t ) ±) cxp(i£ 
Hj o {kR)Ht(nkR) 



— <j)± — 7r/2)), and 
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Here, to = kR sin (9 — <fi±) where the <j>± are the two 
distinct points on the boundary where Snell's law is sat- 
isfied for the incident and reflected ray (see Fig. la, only 
incident ray is shown). 
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Thus Eq. (||) has a simple physical interpretation: 
{a m } describe the relative weight for each angular mo- 
mentum component in the closed resonator, while the 
e m are refraction amplitudes describing the probability 
amplitude for refractive escape from angular momentum 
m in the direction 8. If the interference terms are ne- 
glected in Eq. (p|) t he result is essentially the ray model 
of references ]9|-|12[ generalized for an arbitrary initial 
state, and provides a more rigorous justification for that 
model. However if one maintains (as we do in Fig. 1) 
the generality of Eq. (^|) including the cross-terms then 
the interference effects neglected in the ray model are 
captured. 
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FIG. 1. Counterclockwise from top right: (a) Schematic of 
the classical ray emission, (b-d) Far-field emission intensity 
for a low angular momentum ("bouncing ball") mode with 
nkR ~ 150., refractive index n = 1.5 (b), 2.5 (c), 3.3 (d) and 
deformation of e = 0.106. The black, red and green curves 
show the exact, perturbative (j^ and semiclassical (stationary 
phase) ([]) emission intensities respectively. 

In Fig. 1 we show a test of the emission directionality 
of Eq. (|8|), comparing to an exact numerically-generated 
resonance in the resonator with 10% (e s» 0.1) quadrupo- 
lar deformation R(<fi) = Rq(1 + ecos(2</>)), for different 
values of the refractive index. The exact result (black 
curve) is very well reproduced by the simple perturba- 
tive result of Eq. (||) (red curve) over the whole range of 
indices of refraction, whereas the semiclassical approxi- 
mation (Eq. (^) , the green curve) works well for the lower 
indices but not so well for larger index. The origin of this 
trend is the neglect of the "non-classical" escape pro- 
cesses mentioned above, which become more important 
at high index. These tunneling processes are similar to 
the "ghost" contributions in periodic orbit theory |l6| ] 
and can be included by going beyond the real saddles 
in the stationary-phase approximation as is done in that 
context. We defer this refinement to future work. 

Our theory reduces the calculation of T and 1(9) to 
solving for the eigenstates of the appropriate closed res- 
onator and substituting into Eqs. (H) and (^). In many 



cases it is possible either to evaluate those eigenstates an- 
alytically (e.g. for island or torus states) or to model their 
statistical properties and hence obtain information about 
lifetimes and emission patterns in chaotic resonators. As 
an example, we shall consider the eigenstates localized at 
the stable islands of the mixed phase space of an ARC. 
Such eigenstates are of key importance for the recently 



studied deformed semiconductor microlasers 13 1, where a 
factor of 10 3 enhancement of lasing power was attributed 
to the emission from a mode localized at the stable bow- 
tie periodic orbit (see inset in Fig. 2a). 

We obtain the wave function at the boundary of the 
resonator using the semiclassical theory developed in 
Ref. fl7|| . The quantization on stable islands yields har- 
monic oscillator wave functions; the corresponding oscil- 
lator frequency is determined by the classical dynamics 
in the vicinity of the periodic orbits. For the ground state 
resonance, localized at the center of a chain of stable is- 
land, we obtain the width 



r 



Ro 



20rn 2 *-f R,j, cos 2 x M 



dxe^ 



x — n sin 



(9) 



where Xu is t ne periodic orbit incidence angle and R^ 
the distance from the bounce point fj,. The parameter 
b measures the "spread" of the wave function in angu- 
lar momentum. Choosing as phase space coordinates the 
arc lenth along the boundary and the transverse momen- 
tum sinx, one finds b = (ka^)/ \J 1 + aj t (r) — /) 2 , where 

/ = SpR''/R fi + (R'^/Rfi) 2 / s^, the angular momentum 
I corresponds to the refracted ray, related to the periodic 

orbit [i by sinx<> = «sinx M , and s M = \j (kR^) 2 — P. 

The parameters a M = (4 — Tr [m^-]) 1 / 4 (2rn,2i) 1 ^ 2 and 
r}^ = {1TI22 — Tin ) / m 2i arc defined in term of the mon- 
odromy matrix Q m^. We note that our semiclassical 
method is justified as long as the area in the phase space 
"covered" by the eigenstate, is significantly smaller than 
the size of the islands. 

Similar to the calculation of T, we derived a closed 
semiclassical expression for the emission intensity 1(0). 
This calculation will be presented elsewhere [Q , one re- 
sult for the case of the bow-tie resonance is shown in Fig. 
2b. 

The comparison of the semiclassical emission inten- 
sity and width with the exact calculation is presented in 
Fig. 2. Note, that, as indicated by the Husimi projection 
of the exact eigenstate (panel (a)), the phase space area 
of the island is of the order of the effective fi = 1/nkRo- 
The resulting leakage of amplitude from the island should 
therefore lead to deviations of the actual intensity from 
the semiclassical prediction. However, the general struc- 
ture of the far-field emission pattern (panel (b)) is well 
reproduced by the semiclassical result. 
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FIG. 2. The bow-tie resonance, its emission pattern, and 
its width, (a) Husimi projection (color scale - red denotes 
high intensity) of bow-tie mode at n = 1.5, nkR = 102., 
for quadrupolar defomation e = 0.168 in (</>, sinx) coordi- 
nates, superimposed on the corresponding Poincare Surface 
of Section. The white dashed line indicates the threshold for 
classical emission, sinxc = i/n, above which only tunneling 
escape occurs. Inset: Schematic of bow-tie periodic orbit, (b) 
Far-field emission: Exact numerical results (red) , and the an- 
alytic semiclassical result for this island state (green). The 
formula leading to the green curve is not given in the text 
due to space limitations, it will be derived elsewhere Jl8| . 
(c) Comparison of the exact width of the bow-tie resonance 
(circles and diamonds) with the semiclassical (green line) ex- 
pression (^|) and classical (blue line) prediction. Black circles 
and red diamond correspond to bow-tie resonances close to 
nkRo = 45. and nkRo = 102. respectively. 

In panel (c) we compare the results for the resonance 
width. In addition to the exact data points (circles and 
diamonds) and the semiclassical prediction of Eq. @ (red 
curve), we also show the classical result T = (2Rq/L) lnr 
(blue curve), based on the Fresnel reflection coefficients r, 
calculated at the incidence angle of the periodic orbit (L 
here is the total length of the periodic orbit). The simple 
classical model works very well at index near unity since 
it is not perturbative in T/A which is not very small in 
this case; whereas the semiclassical theory of Eq. ([)]) is 
perturbative and shows a small but visible discrepancy in 
this range. However once the index becomes large enough 
that a ray on the bow-tie orbit will be totally-internally 
reflected the classical model gives the unphysical result 
r = 0, whereas our semiclassical model continues to de- 
crease smoothly and in good agreement with the numer- 
ical data. 
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In summary, we have developed a theory of resonance 
lifetime and emission intensity in nonintegrable dielec- 
tric resonators. The theory is in a good agreement with 
numerical data, has a simple physical interpretation in 
terms of refractive emission, and gives non-trivial predic- 
tions for the lifetimes and emission patterns in asymmet- 
ric resonant cavities. 
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